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INTRODUCTION 


In  recent  years  much  research  has  been  devoted  to  developing  realistic 
constitutive  equations  to  describe  complex  material  behavior  such  as  cyclic 
plasticity  (refs  1-9).  The  majority  of  these  works  are  along  the  lines  of  the 
classical  theory  of  plasticity.  However,  some  attempts  have  also  been  made  to 
establish  new  and  independent  theories.  The  endochronic  theory  developed  by 
Valanis  (ref  8)  is  based  on  the  notion  of  intrinsic  time  and  thermodynamic 
theory  of  internal  variables.  The  original  definition  of  intrinsic  time  has 
led  to  difficulties  in  cases  where  the  history  of  deformation  involved 
unloading.  Valanis  (ref  9)  has  since  introduced  a  new  concept  of  intrinsic 
time  to  overcome  these  difficulties.  The  new  theory  has  been  successfully 
applied  to  describe  the  cyclic  hardening  phenomenon  under  uniaxial  loading 
(ref  10)  . 

In  this  report  the  modified  version  of  the  endochronic  theory  of 
plasticity  is  applied  to  the  problem  of  a  cylindrical  bar  subjected  to  cyclic 
fully-reversed  torsional  loading.  The  governing  equations  of  integral  form 
are  presented  for  pure  shear  def ormation.  Analytical  techniques  are  employed 
in  the  solution  of  these  equations.  The  stress  distributions  in  a  cylindrical 
bar  at  different  stages  of  loading,  unloading,  and  reloading  are  calculated. 
The  solution  for  this  problem  based  on  a  numerical  iterative  technique  was 
reported  recently  (ref  11).  A  different  approach  and  additional  numerical 
results  are  presented  here. 


References  are  listed  at  the  end  of  this  report. 
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ENDOCHRONIC  THEORY 


According  to  the  modified  version  of  the  endochronic  theory  of  plasticity 
(ref  9) ,  the  intrinsic  time  £  for  the  one-dimensional  case  is  defined  by 


dt 


dc  =  |d Y  -  M  —I 


(1) 


where  is  a  positive  scalar  such  that  0  <  <  l,  t  and  y  are  shear  stress 

and  strain  respectively,  and  is  the  shear  modulus.  When  k1  =  0,  the 
original  form  of  intrinsic  time  is  recovered;  but  when  k^  =  1,  flay-  (t/p0) 
is  the  plastic  strain,  and 

dc  =  |dtt|  (2) 

Equation  (2)  is  used  throughout  this  investigation  so  that  the  concept  of 
plastic  strain  may  be  incorporated  into  the  present  development. 

In  the  case  of  pure  shear,  the  governing  constitutive  equation  of 
integral  form  is  given  by 


rz  d  Q 

=  1J0  1  p(z“Zf )  dz f 

o  dz’ 


where 


(3) 

(4) 


u0p(z)  =  y0p0  6(z)  +  Pie“az  +  P2 
in  which  p0,  a,  p  ^ ,  and  P2  are  material  parameters;  6(z)  is  the  delta 
function. 

A  more  general  form  of  the  intrinsic  time  measure  involving  the  strain 
rate  effect  is  proposed  in  Reference  10  as 

dc  =  k(  |a|)  jda|  (5) 

where  k,  the  strain  rate  sensitivity  function,  is  a  function  of  Q.  In  this 
case,  Eq.  (3)  can  be  written  as 

,z  1  dC 

(6) 


z  1  d  c 

=  UQJ  p(z-z') [±  -]  -  dz' 

o  k  dz 1 
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The  constitutive  Eq.  (6)  is  applied  here  to  describe  the  material  behavior  at 
constant  plastic  strain  rate.  The  intrinsic  time  z  is  related  to  £  by  the 
following  time  scale 

dC 

(7) 


dC 

—  =  f(C) 
dz 


where  f(£)  describes  isotropic  hardening  and  is,  therefore,  termed  the 
hardening  function.  In  this  report,  the  form 

f(z)  =  c  -  (c-l)e~^z  (8) 

is  used  because  of  its  simplicity  and  its  proved  usefulness  in  cases  of  cyclic 
loading  (ref  10).  The  parameters  c  and  3  are  material  constants. 

CYCLIC  SHEAR  RESPONSE 

By  using  the  hardening  function  given  in  Eq.  (8),  the  constitutive 
equation  for  loading  with  constant  plastic  strain  rate  is  written  as 

T  =  (M0/k)/  p(z-z' )f(z')dz'  (9) 

o 

where  p(z)  is  defined  by  Eq.  (4)  and  k  is  a  constant.  Integrating  Eq.  (9), 
the  following  explicit  result  is  obtained: 

T  =  (Ty/k)f(z)  +  (ui/k)(g(z)-g(o))  +  (y2/k)h(z) 


where 

Ty  =  the  yield  stress 

h(z)  =  cz  +  3-1(c-l) (e~^z-l) 

§(z)>  g(o)  are  the  values  of  the  function  g(z')  evaluated  at  z* 
respectively,  and 

g(z')  =  (c/a)e-a(z-z' )  -  (c-l)/(a-3)e_az+(a-3)z' 


(10) 

(ID 

(12) 


=  z,  o, 


(13) 
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Equation  (10)  is  now  the  response  function  for  loading.  If  unloading  occurs 
when  the  plastic  strain  reaches  Oj  (or  z  =  z{)  ,  then  the  response  function  for 
z  >  can  he  obtained  as 

t  =  -(xy/k)f(z)  +  (M1/k)[-g(z)  +  2g(Zl)  -  g ( o ) ]  +  (p2/k)[-h(z)  +  2h( zx) ]  (14) 

If  Eqs.  (10)  and  (14)  are  examined  at  z  =  z±~  and  zf1",  respectively,  a  drop  in 
stress  of  2(xy/k)f(z^)  results  during  elastic  unloading. 

If  reloading  takes  place  after  unloading  and ,  assuming  that  reloading 
occurs  at  £2  =  S22  (°r  z  =  z2>  where  z2  >  zj),  the  response  function  for  z  >  z2 
is  obtained  as 

T  =  (xy/k)f(z)  +  (p1/k)[g(z)  -  2g( z2)  4-  2g ( z L )  -  g(o) ] 

+  (M2/k)[h(z)  -  2h( z2)  +  2h(z1)]  (15) 

When  Eqs.  (14)  and  (15)  are  examined  at  z  =  z2~  and  z2+,  respectively,  a  jump 
in  stress  of  2(xy/k)f(z2)  is  again  obtained  during  elastic  reloading. 

This  procedure  of  obtaining  a  theoretical  expression  for  each  part  of  the 
cyclic  loading  process  may  be  continued.  Thus,  a  general  expression  of  the 
response  function  for  the  constant  total  strain  amplitude  cyclic  torsion  test 
may  be  found  to  be 

x  =  (-l)N(xy/k)f(z)  +  (u1/k)G(z)  +  (u2/k)n(z)  (16) 

where 

N 

G(z)  =  (“l)N[g(z)-g(zN) ]  +  \  (-l)n+1[g(zn)  -g(zn_1)]  (17) 

n=  1 

N 

>Kz)  =  (-l)N[h(z)  -  h(zN)  +  l  (-l)n+1[h(zn)  -  h(zn_L)  (18) 

n=l 

Note  that  N  ^  0  is  the  number  of  half —cycles ,  odd  for  unloading  and  even  for 
reloading . 
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After  many  cycles  when  the  values  of  z  become  large,  the  hysteresis  loop 
of  the  stress-strain  curve  will  approach  a  steady  state.  If  Eq.  (16)  is 
examined  at  z  =  z^“  and  z  =  ZN*">  a  drop  or  jump  in  stress  of  magnitude 
2(Ty/k)f(zN)  results,  which  corresponds  to  the  elastic  response  upon  the 
reversal  of  the  loading  or  unloading  direction.  When  z^  is  sufficiently 
large,  the  jump  or  drop  of  stress  becomes  a  constant  value  2c  Ty/k  at  the 
steady  state. 


CYLINDRICAL  BAR  UNDER  TORSION 

For  a  cylindrical  bar  under  torsional  loading,  the  external  torque  is 

ra 

Ts  =  2tt  /  tr2dr  (19) 

o 


where  r  is  the  current  shear  stress  corresponding  to  location  r,  and  ra  is  the 
radius  of  cross-section.  Geometric  considerations  show  that  radial  lines  have 
to  remain  straight  after  deformation.  Thus,  one  concludes  that 


Y  =  (Ya/ra)r 


(20) 


where  Ya  is  the  strain  at  the  outermost  fiber.  Since  a  yield  stress  is 
introduced  in  Eq.  (4),  an  elastic  core  always  exists  during  deformation  whose 
radius  re  is  given  by 


Ty  ra 


(21) 


^o  Ya 

If  the  experiment  is  strain-controlled  with  strain  at  ra  varying  between 
~Ya  and  +ya,  and  with  Ya  in  the  plastic  range,  then  the  torque  can  be  computed 


as 


TS  =  -  Tv  re3  +  2ir(ra/ Ya) 3/  TY2dy 


(22) 
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where 


Y e  =  (Ya/ra)re  ,  Y  =  &  +  t/u0  (23) 

t  and  are  given  by  Eqs.  (16)  and  (18),  respectively. 

Now  that  y  is  known,  the  values  of  z  and  t  at  each  fiber  can  be 
calculated  and  used  in  Eq.  (22).  Note  that  there  exists  an  explicit 
expression  for  all  cycles 

1 

dy/dz  =  Q'( z)  +  —  x'(z)  =  p(z)  (24) 

where 

fi'(z)  =  (-1)N  f(z)  ,  f'(z)  =  B(c-l)e“Pz 

t'(z)  =  (-l)N[ty/k)f ' (z)  +  (P2/k)f(z)l  +  ( Ml/k) T* ( z) 

N 

T'(z)  =  (-1)N[F'zz  -  F'Nz]  -  l  (-l)n[F'nz  -  F'(n-i)z] 

n=l 

F'zz  =  [8(c-l)/(a-|3)]e“^z  ,  F'nz  =  -aFnz 

and 

F’nz  =  "  e“a<z”zn)  -  (--i)e"az  +  (a"P)zn  (25) 

a  a- 3 

The  integral  in  Eq.  (22)  can  be  replaced  by 

Y  a  za 

/  TY^dy  -  /  Ty2p(z)dz  (26) 

^e 

where  za  can  be  calculated  by  Eq.  (23)  with  known  value  Ya.  Now  the  numerical 
integration  becomes  very  easy  because  the  values  of  the  integrand  can  be 
evaluated  directly. 
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NUMERICAL  RESULTS  AND  DISCUSSION 


To  apply  the  developed  model,  the  material  constants  (a,  3,  c,  Ty,  pQ , 
p  i ,  P2)  in  the  theory  have  to  be  determined.  These  material  constants  can  be 
determined  if  the  cyclic  shear  stress-strain  curve  for  the  material  has  been 
obtained  experimentally .  The  usual  procedure  is  to  perform  a  cyclic  torsion 
test  using  a  thin-walled  tubular  specimen.  Such  a  test  for  the  annealed  AISI 
4142  steel  was  carried  out  by  the  Plasticity  Research  Laboratory  at  the 
University  of  Iowa.  The  values  of  constants  were  then  used  to  predict  the 
results  for  a  solid  bar  test.  The  theoretical  and  experimental  results  were 
in  reasonable  agreement  (ref  1L).  This  real  material  did  not  show  any 
appreciable  amount  of  cyclic  hardening. 

For  purpose  of  investigating  the  implications  of  the  developed  model,  a 
hypothetical  material  with  appreciable  cyclic  hardening  behavior  was  studied. 
The  shear  stress-strain  behavior  of  such  material  under  fully-reversed 
torsional  loading  is  presented  in  Figure  1.  The  material  constants  were 
determined  as:  a  =  1000,  3  -  50,  c  ~  1.5,  iy  =  104  psi,  p0  =  107  psi,  p^  = 
4xl06  psi,  p?  =  0.  A  steady  loop  was  established  after  a  few  cycles.  The 
same  set  of  constants  was  then  used  to  predict  the  stress  distribution  in  a 
cylindrical  bar  subjected  to  cyclic  torsional  loading.  We  carried  out  the 
computational  process  seven  cycles  after  initial  loading.  The  numerical 
results  are  presented  here  for  the  initial  loading  half-cycle  and  two  cycles 
of  unloading  and  reloading. 

Figure  2  presents  the  numerical  results  for  the  torque  as  related  to  the 
shear  strain  at  the  outermost  fiber.  The  distribution  of  stress  in  the 
cross-section  at  different  magnitudes  of  torque  during  the  initial  loading 
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half-cycle  is  presented  in  Figure  3.  The  corresponding  shear  strains  at  the 
outer  fiber  are  0.1,  0.3,  0.6,  and  1.0  percent,  respectively.  Notice  that  the 
outer  fiber  is  the  first  one  to  yield  at  Ya  =  0.1  percent;  subsequently  as 
more  torque  is  applied,  the  radius  of  the  elastic  inner  core  gets  smaller. 

Also  notice  that  the  rate  of  hardening  for  each  fiber  decreases  which  is,  of 
course,  in  accordance  with  strain  hardening  phenomena. 

Figure  4  presents  the  distribution  of  the  shear  stress  in  the  bar  at 
three  stages  of  the  first  unloading  half-cycle.  The  top  and  bottom  curves 
correspond  to  the  beginning  and  end  of  the  unloading  stages,  i.e. ,  Ya  -  -  1 
percent.  The  corresponding  values  of  torque  are  1256,  -1387  lb-in., 
respectively.  The  middle  curve  represents  the  residual  stress  distribution 
when  the  applied  torque  is  equal  to  zero.  A  few  iterations  are  needed  to 
reach  this  state  and  the  residual  shear  strain  at  the  outermost  fiber  is 
0.378  percent.  Figure  5  presents  the  distribution  of  the  shear  stress  in  the 
bar  at  three  stages  of  the  first  reloading  half-cycle,  i.e.,  Tg  =  -1387,  0, 
1461  lb-in.  The  middle  curve  represents  the  distribution  of  residual  stress 
when  the  torque  is  equal  to  zero.  Reverse  yielding  occurs  within  the  outer  12 
percent  of  the  section  when  Tg  =  0  during  unloading  and  reloading.  The 
distribution  of- residual  stress  during  unloading  is  quite  different  from  that 
during  reloading  as  shown  in  Figures  4  and  5  for  the  first  cycle.  Similar 
results  for  the  residual  stress  distribution  during  the  second  unloading¬ 
reloading  cycle  are  shown  in  Figure  6.  The  solid  curve  is  the  favorable  one 
if  the  applied  torque  during  service  is  in  the  same  direction  as  the  initial 
loading . 
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Figure  2,  Torque  vs.  Outer  Shear-Strain  Curve. 
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radius  r/ra 

Figure  3,  Stress  Distribution  During  Initial  Loading* 
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radius  r/ra 

Figure  4.  Stress  Distribution  During  First  Unloading. 
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Figure  5.  Stress  Distribution  During  First  Reloading 
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Figure  6,  Distribution  of  Residual  Stresses  in  Second  Cycle. 
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